Addition Theorems Via Continued Fractions 

Mourad E.H. Ismail Jiang Zeng 

Department of Mathematics Universite de Lyon, Universite Lyon 1 

University of Central Florida Institute Camille Jordan, UMR 5028 du CNRS 
Orlando, FL 32816 69622 Villeurbanne 

USA France 

February 2, 2008 



Abstract 

We show connections between a special type of addition formulas and a theorem of 
Sticltjes and Rogers. We use different techniques to derive the desirable addition formu- 
las. We apply our approach to derive special addition theorems for Bessel functions and 
confluent hypergeometric functions. We also derive several additions theorems for basic 
hypergeometric functions. Applications to the evaluation of Hankel determinants are also 
given . 

Mathematics Subject Classification Primary 33D15, 33 C15, Secondary 30E05, 05A15. Key- 
words Addition theorems, orthogonal polynomials, continued J-fractions, g-orthogonal poly- 
nomials, Askey- Wilson polynomials, Bessel and confluent hypergeometric functions. 

1 Introduction 

An algebraic addition theorem for a function / is an identity of the form 
(1-1) P(f(x),f(y),f(x + y))=0 

for some polynomial P in three variables. Weierstrass proved that an analytic function satis- 
fying an algebraic addition theorem is a rational function in z, a rational function in e^ z for 
some A, or an elliptic function, Chapter 13]. This notion is too restricted to be useful in 
the theory of special functions. In general a family, say (f)\ , of special functions satisfies an 
addition formula if there is an elementary continuous function A of three variables x, y, t and 
an expansion in terms of a family of special functions Vv such that the expansion coefficients 
factor as products in x and y. In other word we have 

(1.2) <f> x (A(x,y,t)) =^C(A, / u)^(x)<(y)^(t), C(A,/x)eC. 
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Recall the definition of a Bessel function 

(1.3) - m!r(m + z/ + 1) J - r>TI) ° Fl V 1/ + 1 —) 

and the modified Bessel function 

oo 

(1.4) J„(x) := / 2 J,(e^ 2 x) = £ 



(x/2f+ 2m 
m! T(m + zv + 1) 

m=0 



One important addition theorem is the addition theorem for Bessel functions, 

for i/=£ 0,-1,-2,- ••, where w := (z 2 + Z 2 -2zZcos</>) 1/2 , [7, (7.15.30)], [26J. The polynomials 
{C^(x)} are the ultraspherical polynomials. The special case (ft = it is 

h a Ju{x + y) Tjy) A, , , (-l) w (2i/) w 

(L6) = (^F ^ ( + n) — ^ — " +n(a;) ^ +n(y) ' 

since C^(— 1) = (— l) n (2zv) n /n!. The addition theorems we will encounter in this work are of 
the type (fL6|) . 

This work arose from an attempt to understand the Stieltjes- Rogers theorem of continued 
J- fractions, see Theorem 12.11 It is clear that (|2.12p of Theorem 12.11 is an addition theorem of 
the type (fTBjh 

We decided to explore ^-analogues of Rogers' addition formula (|2,12p of Theorem 12.11 and 
to compute the functions Qj(x),j = 0, 1, . . . for specific continued fractions, since the theory 
of orthogonal polynomials, especially the recently discovered one, provide a rich source of 
continued J- fractions. We discovered two g-analogues of Theorem 12.11 They are Theorems 14.11 
andlPl 

In this work we establish additions theorems for many special functions. To the best of 
our knowledge only (|1.6|) and (|3.8|) are known. We offer three different techniques of proof and 
provide at least one example of each technique as illustrations. We realize that it is possible 
to use fewer techniques to achieve the same goals but we believe there is merit in utilizing as 
many different ideas as possible. One approach uses the plane wave expansion, [151 (4-8.2)], 



(17) „_y- r(q + /3 + n+l) / (3 + n+l 

71=0 

for a > —1, (3 > —1, and its special case [El (4.8.3)], 

00 

(1.8) e^ = r(z,)(y/2)-^(z, + n)/, +n (y)C(x), v > -1/2. 

n=0 



2y) Pt>®(x), 
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The polynomials {P^ 1 '^ (x)} and {C^(x)} are Jacobi and ultraspherical polynomials, respec- 
tively. The expansions (|1.7|) ~ (|1.8|) are instances of the Fields and Wimp expansions [9], see 
also [8], [24J. Other techniques use Rodrigues type formulas followed by integration by parts, 
and connection coefficient formulas. 

We shall follow the notation and terminology in [I], |15j . and |12j . In particular we shall 
use the Rogers connection coefficients formula for the continuous g-ultraspherical polynomials 
{C n (x;0\q)}, 

(1.9) C B (x; 7 | ? )=g (?)fc(A _ fe r3 ^C n „ 2fc( x; /?!,), 

|15l p. 330], and the facts that 

(1.10) CT n (x) = C7i(o;) = C7 n (x;g| 9 ). 



2 Preliminaries 

Given a moment sequence {fJ- n }, we define the linear functional C : x n i— > /i n on the vector space 
of polynomials C[x]. We shall always assume //o = £(1) = 1- Then the monic polynomials 
P n (x) orthogonal with respect to the C or the moment fi n satisfy the following three term 
recurrence relation (the spectral theorem for orthogonal polynomials [15, Chapter 2]): 

(2.1) P n +i(x) = (x-b n )P n (x) -A„P„_i(x), n>0, 

where \qP-i(x) = and Po(x) = 1. We shall always require the functional to be regular, [5], 
which is equivalent to demanding that A n 7^ for all n,n > 0. The orthogonality relation is 

(2.2) C(P m P n ) = X1X2 ■ ■ ■ X n ^m,n- 



The moment sequence is related to the coefficients b n and A n by the following identity: 



Define the determinants 



A 



i,n 



Mo 

m 



Mi 

M2 



Mi-1 Mi 
Mn Mn+1 



Mi 
Mi+1 

M2i-1 
Mfi+i 



£>n(x) 



Mo Mi 
Mi M2 



Mn— 1 Mn 
1 X 



Mn 
Mn+1 

M2n-1 



In particular, let 
(2.4) 



A, = A r 



Xr, 



A 



n,n+l • 



Then P n {x) = (D n —\) D n (x) is the monic orthogonal polynomial sequence for C. 
It is easy to see that 



(2.5) 
(2.6) 

Therefore 
(2.7) 
and 
(2.8) 



C{x n P n {x)) 
C(x n+1 P n (x)) 

\„ 

b n = 



D„ 



D n -1 
Xn 

D n -i 



A-nAn-l • • • Al, 

Ai(6 



+ -" + bn). 



£[P n 2 (x))] D n _ 2 D r , 



C[xPl{x))\ _ Xn__ Xn-1 
D n -1 



C[P^(x))] D h 
The next theorem is the backbone of this work. 

Theorem 2.1. Define the Stieltjes tableau of entries Hi^ n (i,n > 0) by 



(2.9) 



Hi,n = for i < and i > n; 
#n n = 1 / or a ^ re > 0; 

-ffj.n = Hi-l,n-l + biPli^n-x + Aj+l-E/j+i n _i< 



T/ien £/ie generating function Y2n>o Ho,n xU has the continued fraction expansion f|2.3|) i/ and 
only if, for any two nonnegative integers k,£ >0, the following convolution identities 



(2.10) 



H, 



o,k+e — Hq^Hqi + XiHij,Hi i H h Ai • • • XjHj k Hj £ + 
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hold. Moreover with the exponential generating functions of {Qj(t)} 

00 t n 

(2.11) Qj{t)=Ys H i,n- v 

n=j 

the convolution identity (|2.10p is equivalent to the addition formula 

00 

(2.12) Q (x + y) = J2 Ai • • • A n Q n (x)Q n (y). 

n=0 

Wall [23] points out that the first part of Theorem 12.11 is due to Stieltjes but the addition 
theorem (|2.12p is due to Rogers. For a proof and references see [231 Section 53]. 
It is important to note that the Hj^'s are the connection coefficients in 

n 

(2.13) x n = Y J H hn P 3 (x). 

Observe that the addition formula fl2 . 12f) is equivalent to 

fotA\ xh (x) - yh (y) ^ 

(2.14) = > Ai • • • A„ h n {x)h n (y), 

x — y ' 

where 

00 

n=j 

In general (I2.10p implies 

00 m 00 

(2.15) H 0>™ C^dm-jt™^ = Y, A l • • • X n Qn(s)R n (t), 
m=0 j=0 n=0 

where 

00 00 

(2.16) Q n {s) = I'':,'-,-'- Rn{t) = ^ HnjdjV . 

j=n j=n 

We may define a generalized translation operator (GT) on polynomials by 

m 

(2.17) (GT) s x m = c jS id m ^t m -i. 
One can extend (GT) S to formal power series by linearity. 



Remark 2.2. According to the Flajolet-Viennot theory [10\ 125] we can interpret Hi^ n in the 
Stieltjes ' tableau as follows. Let us attach weights to the steps of a lattice path at level i (i > 0) 
of a Motzkin path in the following way: 

w(/) = 1, w(-) = hi and w(\) = Aj. 

Let To^i(n) be the set of Motzkin paths from level to level i of length n. Then we have the 
following interpretation: 

AiA 2 • • • \H i>n = 10(7). 

76r^ (n) 

This provides a combinatorial interpretation of ([2, 101) . 

Let {Pj(x)} be the monic orthogonal polynomials with respect to the moment sequence 
i?0,n an d £ be the functional C{x n ) = H n . Then it follows from f)2 . 13|) that 



Note that hi = H i+ i, i+2 - H i>i+ i and 

(2-19) //,.■ t^V- 

Therefore the addition formula (I2.12D generalizes the Hankel determinants in (12, 4p . We refer 
the readers to [3j [13j EH EH [28] for the application of orthogonal polynomials to the 
computation of Hankel determinants. 

The following theorem gives another interpretation of the function Qo(x) for which our 
technique will derive an addition theorem. 

Theorem 2.3. Assume that {P n (x)} are orthogonal with respect to a positive measure /i with 
compact support contained in {z : \z\ < r}. Then 

(2.20) AiA 2 . . . XjQj(t) = <b e t/z F n (l/z 



z 



where 



f P (u) 

(2.21) F n (z) := / ^^dfi(u), z $ supp{//}. 
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Proof. Let /x n = j R x n dn(x). The right-hand side of (|2.20p is 

I e t /z f PM_ d (u) v ^ fwp^u) d M (t 



n=Q k=0 



Jr^ n\ J R 



and the theorem follows from (12. 18H . 

The function F n {z) is related to the function of the second kind [15]. 
Throughout this work we will use the Heine transformation 

(2.22) 2 MAB,C,Z) = ^L^l^ 2M c/B,Z;AZ;q,B) 
[12\ (III.l)] and the 2 4 > i to 202 transformation 

(2.23) 2 MAB;C;q,Z) = 5 f 00 2 fc(A C/B; C, AZ; q, BZ). 

m (ni.4)]. 



3 Ultraspherical and Jacobi Polynomials 

The ultraspherical (or Gegenbauer) polynomials are 

(2t/) w / -n, n + 2v l-x 
° n{x) ~ ^T 2jFi V I/+1/2 



The normalized weight function is 



The corresponding orthogonality functional C is defined by 



£(/) = y f(x)w(x)dx. 



The monic ultraspherical polynomials {-P n (x)} and the A's are 



4(i/ + j-l)(i/ + j)" 
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Moreover 

„ / s „ / n n(n + 2v — 1) „ 

X P n {x) = P n +l(x) + 7r -^ ' / n-l4 

Therefore (|2.18|) when f > — 1/2 implies 



A\A2 ■ ■ ■ Aj 

^ n -r(,)(t/2)-£ ^f> + n)U(i)C:(,)P 4 (x)j 

= r(i/)(t/2)-"?^(i/ + i)7 I ^(t). 



Therefore 



(3-D Qi(t) = 2T ^)^ +1) ^(t), 

It is straightforward to see that addition formula f)2 . 12|) in the present example is equivalent 
to CPU) . 

Next we consider Jacobi polynomials. The normalized weight function is 
(3.2) w(x) = (l-x) a (l + xfA(a,P), A(a,0)~ I <o + ./+>) 



2«+/J+ir(a + l)r(/3 + 1) ' 
In this case the functional to be considered C is 

£(/) = y f(x)w(x)dx. 

The monic Jacobi polynomials {P n (x)} are defined through 

p(^)(x) = (n + a + / ^ + 1) " P ra (^ 
n v i 2 n n! 

so that 

4n(n + a)(n + /3)(n + a + /3) 



(3.3) A n 



(2n + q + /3 - l)(2n + a + /3) 2 (2n + a + /3 + 1) ' 
As in the case of ultraspherical polynomials we can use (jl.7p and the result is 



(3.4) Qi(x) = — — 1F1 



a + (3 + 2* + 2 



2i 
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As an example of the use of Rodrigues formulas we give another derivation of (|3.4p . The 
Rodrigues formula for Jacobi polynomials is \15\ (4.2.8)] 



(3.5) (l-x) a (l + x)Ppt^(x) 



(-l) n d n 
2 n n\ dx™ 



n+/3 



Therefore 



Qj(t) 



Vj\A(a,(3) 



X 1 ...X j (a + (3 + j + J_ x 
(-iyA(a,P) 



\l...\j(a + P + j + l)j J_i dxi 



xt 



(l-x) n+a (l + x) 



(l-x) a (l + xfP^(x) dx 

d? 



(l-x) j+a (l+x) j+/3 



dx 



t j A(a, (3) 



(l-x) j+a (l + x) j+0 dx, 



\ l ...\ j (a + (3 + j + J_ x 
after integration by parts. Taking into account the integral representation (6J 6.5.1)] 

r(c) 



(3.6) iFi(a;c;z) 



[ e zu u a - l {l-u) c ' a ~ l du, Rec>Rea>0, 
J o 



T(a)T{c-a) 

(|3.3p . and (|3.2[) . we see that the last expression for Qj(x) reduces to 



(3.7) 



Qj(t) = %e t 1 F i 



a + j + 1 
a + (3 + 2j + 2 



21 



The equivalence of the representations (|3.4p and (|3.7p follows from the transformation [15(. 
(1.4.11)] ' 



1*1 



c — a 
c 



= e\F l 

This analysis establishes the following theorem, which is the main result of this section. 
Theorem 3.1. We have the addition theorem for the confluent hypergeometric functions 



1*1 



(3.8) 



a + 1 

a + P + 2 

Xi.Fi 



/ ! | £ 

n=0 

a + n + 1 



(a + l) n (f3 + l) n (a + f3 + l) n (tsf 
(a + (3 + l) 2 „(a + /S + 2) 2n n! 



* I 1*1 



a + n + 1 
a + P + 2n + 2 



a + j3 + 2n + 2 
When a = /3 = z^ — 1/2 Theorem 13.11 reduces to (jl.6p since 
(3.9) e~ x iFi(i/ + 1/2;2i/ + 1;2x) = 1> + l)(2/x) l/ / 1 ,(x), 
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[6J (6.9.10)]. Moreover both ([33|) and ([377D also reduce to ([3dD . 

Note that Theorem 13.11 and (|3.1|) can be proved from Theorem 12.31 and the facts 




(z -!)"+! 2 1 Va + /3 + 2n + 2 1-zJ 



see (4.4.1) and (4.4.6) in [15j . Moreover the integral 



/_ 



i 



e xt {l-x) a {l + xf P^\x)d. 



x 



i 



can be evaluated from the plane wave expansion (jl.7jl instead of the use of Rodrigues' formula. 

It must be noted that (|3.8p coincides with formula (42) of Burchanl and Chaundy [I]. It is 
also a limiting case of formula (50) in [3]. The latter is stated in [6], see the first unnumbered 
formula after (7) in Section 2.5.2. Indeed if we replace z by z/b and £ by (,/b and let b — > oo, the 
above mentioned formula reduces to our (|3.8|) . The terminating case a = —m — 1, (3 = j + m of 
(|3.8p is the case r = of Koornwinder's addition formula for Laguerre polynomials, see (3.3) 
in [19]. Also, this terminating case of (|3.8p is the inverse of 10.12(42) in [7]. 

4 Two g-Addition Formulas 

The (/-binomial formula 




CO 



(4-1) 



yields Euler's g-analogues of exponential formula: 



(4-3) 



(4.2) 




The g-difference operator is defined by 



m-f(qz) 
(1-3)* 
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It is easy to see that 

(4.4) ^{-xt-qU) = 7r A_(-x^;g) c 



(4.5) V q 
Note also that 



(1 - q)i 
1 \ _ P(l - q)-J 



(xt;q) OD J (a;t;g)oo 

Pq~(i) i 



(4.6) V 

(4-7) V'Ji-xt-qU) = 1 *(- x t;q) 00 . 



llq \{xt;q) 00 ) ' (i-?)Ma*r J ;?)c 

1/q ((~ xt > Q)oo) = (it q y 



4.1 First q- Addition Formula 

Define the ^-translation operator T s ^ q by 

(4.8) T yiq x n = {x + y){x + yq)---{x + yq^ 1 ). 
Extend this by linearity on functions f(x) = J2^Lo a nX n , so that 

oo 

T y ,qf{x) = an ^ x + y)( x + yi)---( x + m 11 ' 1 ) 

n=0 

Note that 

(4.9) T y , q -ix n = {x + y){x + y/q) • • • (x + y/q^ 1 ) = %, q y n q 



Define two ^-analogues of Qj(t) by 



t n ~ l± s ^ n q^)t n 

' n't 



(4.10) Qj(t;q) = ^^. n — —, Q^t-q) =J2 H i, 



Theorem 4.1. TTie convolution identity ()2.10p is equivalent to the addition formula 

oo 

(4.11) T Sjq Q (t; q) = ^ A x • • • A n Q„(£; g) Q n (s; g). 

n=0 

Moreover, we have 

(4.12) Qj ^ q ) = ^-^. C ( P / ir) 



Ai---Aj \(a:*;g)c 

(4.13) Q^i; g) = 1 £ (P J (x)(-xt; qU) . 

Ai • • • Aj 
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Proof. The (/-binomial theorem and the definition of T y gives 

n 

T yA x n = (x + y)(x + yq)---(x + yq n ~ l ) = ^ 



fc=0 



q\2) X y 



This establishes the equivalence of (|4.1ip and (|2.10|) . Equations (|4.12p and (|4.13j) follow from 
Euler's formulas K^-KM- □ 



4.2 Little Jacobi Polynomials 

The little g-Jacobi polynomials {p n (x;a, b)} are defined by 

(4.14) p n (x;a,b) = 2 MQ~ n ,abq n+1 ;aq;q,qx). 

[TBI p. 92-93]. The corresponding monic polynomials {Pj(x)} are given by 



(4.15) ^^^ '-"•f^'^'- F.M, 



and 

(4.16) 

Let 



aq 2n-l^ _ gnyj_ _ ggnyj_ _ bq n^(i _ ahq n^ 
(1 - abq 2n - 1 ){l - abq 2n )(l - abq 2n )(l - abq 2n + l ) ' 



(4.17) a = q a , b = 
then the weight function is given by 

w{x;a,f3) = — g— — x , 

and the corresponding orthogonality functional is 

(4.18) £(/) = . =■ ( 2 ° g ' q )°° : f f(x)w(x; a, (3)d q x. 

(a6g 2 ,g;g)oo(l - q) J 

The Rodrigues-type formula for the little g-Jacobi polynomials is 

1 (1 - aVa ja+ (i) , 

<" 9 > »<■*■• « " ^" {W(X] " + + f>) - 
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Therefore 



Qj(t; < 



a iq i (abq,abq 2 ;q) 2j (aq,bq;qU f 1 3)P(x)d x 

[q, aq,bq, abq; q) jiabq^, q; qjooil - q) J 

°—\ w g , ) f 2 ^ aq,b ^ (-l)Jg(a) / (-xt;g) oc u;(x;a,/3)p i (x;a,6)dg2; 

{q,°q;q)j{abq,q;q)oo{i-q) J 

q~ j (abq 2 ;q) 2j {aq,bq;q) 00 (-l) j {l-q) j r 1 



(q, aq,bq; q)j(abq 2 , q; q)^ 
The g-analogue of integration by parts is 



(4.20) 



V q (f(x))w(x)d q x 



(-xt; q) 00 D J 1/q (w(x; a+j,8 + j))d q x. 



f(x)V 1/q (w(x))d q x, 



provided that w(a/q) = w(b/q) = 0, see |T5 S (11.4.9)]. Applying (|4.20p to the last expression 
for Qj(t; q) we obtain 



Qj(t; 



(abq 2 ; q) 2j (aq, bq; q)oot j q^ A (~tq j+n , q n+1 ; g)oo ( a +j+l)n 

/ j I _fl4-r>4-1 4-1 . ~\ ^ 



(q,aq,bq;q) j (abq 2 ,q;q) OQ ^ ( q f3+n+l+j. ? ) 



-tqi,aqi +1 ; q)^q® ( bqi +1 , 



(q;q)j(abq 2+2j ;q)c 



29l 



-tq> 



q,aq 



This shows that 



(4.21) Qj{t- q ) 



t j q® (-tqi,aqi +1 ;, 
lim 



- 3V "'^ (<?;</),■ ^6 (aV+2,, e; ^ 
Now the Heine transformation (|2.22|) leads to 

(-igW +1 ; 



2</>l 



6g^ +1 , e 
-tq j 



q,aq 



j'+i 



"2^1 



6g>' +1 , e 
-tqi 



q,aq 



J+l 



2</>l 



-tg'/e, ag J+1 
a6g 2 J+ 2 



(abq 2 +V,e;q) c 

Therefore the above equation and (|4.2ip establish the basic hypergeometric representation 



(4.22) Qj(t;< 
Similarly we have 

Qj(t;q 



fig (a) 



i<Pi 



ag-?' +1 
26g 2 ^+ 2 



g, -tq- 



(abq 2 ;q) 2j (aq,bq;q) OD P ^ (g 1+ "; g) 00 a"g^+ 1 )" 



{q,aq,bq;q) j (abq 2 ,q;q) 00 ^ (bq n+1+ i ,tq n ; q) c 



t J W +1 ;?)c 



(q;q)j{t,abq 2 i +2 ;q) c 



291 



bqi +1 , t 




g,ag 
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Applying the Heine transformation (|2.22p to the last 291 yields 

t j 



(4.23) 



Qj(t; < 



(?;?). 



-291 



0, aqi +1 
abq 2 ^ 2 



q,t 



Theorem 4.2. For the functional in (|4.18p i/ie functions Qj{t;q) and Qj(t;q) are given by 
(|4.22p and ()4.23p . Moreover we have the following addition formula: 



291 



(4.24) 



aq, —s/t 
abq 2 

*291 



q,t / 

3=0 
0, ag^' +1 
a6g^'+ 2 



(q;q)j(abq,abq 2 ;q) 2j 
,j+i 



q,t i9i 



aq> 
abq 2 ^ 2 



q, -sq J 



Proof. We only need to show that the left-hand side of (|4.24|) is T s ^ q Qo(t, 
the case as can be seen from (14.81). 



This is indeed 
□ 

It is important to note that Theorem 4.2 is a g-analogue of Theorem 3.1. Indeed with 
a = q a , b = q@ and s and t replaced by s(l — q) and t(l — q), respectively, equation (j4.24j) tends 
to <|32H as q -> 1". 

Note that the transformation (|2.23p implies 



(4.25) 



Qj(t;q) 



(q;q)j(t;g)c 



-191 



bqi +1 
abq 2 i +2 



q,q 



i +l at 



Thus the addition theorem (|4.24p has the alternate form 

f(3{aq,bq,abq;q)j (qast) 3 
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(4.26) 



aq, —s/t 
abq 2 



q,t 



3=0 



(q;q)j(abq,abq 2 ;q) 2 j (t;q) c 



Xi9i 



bqi +1 
abq 2 i+ 2 



l,q j+l at\ !9i ( 



aqi +1 
abq 2 ^ 2 



Remark. A different g-analogue of Theorem 3.1 was given by Jackson |171 (55)]. 

4.3 Big g-Jacobi Polynomials 

The monic big g-Jacobi polynomials are 



(4.27) 

Let 

(4.28) 



(abq n+1 ;q) n 
p n (x;a,b,c) = — P n [x). 



W\{x, a, b, c) 



(aq, cq; q) n 

_ (x/a,x/c) c 
(x,bx/c)o. 
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and 

w±(x, a, b, c) (aq,bq,cq,abq/c;q) c 



(4.29) w(x,a,b,c) 



aq(l-q) (q,abq 2 ,c/a,aq/c;q) c 
The corresponding orthogonality functional is 

r-aq 
I cq 

The Rodrigues-type formula for big g-Jacobi polynomials is 



,-aq 

(4.30) L(/)= / f(x)w(x,a,b,c)d q x. 

J cq 



(4.31) Wl (x)p n (x;a,b,c) = } i*L.< D %w 1 (x,aq n ,bq n ,cq n )- 

(aq,cq;q) n q 

Note that 

-acq n+1 (l - q n )(l - aq n ){l - bq n )(l - cq n )(l - abq n )(l - abq n /c) 



(4.32) A r , 



(1 - a6g 2 «- 1 )(l - abq 2n ){\ - abq 2n ){\ - abq 2n + v ) 
So 



(4.33) Ai ■ ■ ■ A 

Therefore 



(— acy q^^ +3 ^ 2 (q, aq, bq, cq, abq, abq/c; q)j 
(abq,abq 2 ;q) 2 j 



Qj(t,a,b,c) = - — — : -h(P j (x)/(xt;q) 00 ) 
Ai • • • \j 



{-iy (1 - qyiqi't 1 ) {aqi +l , bqi +l , cqi +1 ,abqi +1 /c; « 



(q; q)jaq{l - q){q, abq 2 , c/a, aq/c; q) c 

J cq 



x / VqWitx^q^bq^cqi) d ^ 

I cq { x ^i Q)co 

Vq-ljabq 2 ; q) 2j (aq ]+1 , bq^ +1 , cq^\ abq^/c; q)^ f 
a(q;q)j(q, abq 2 , c/a, aq/c; q)oo 31 
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where 



q(l- 

oo 

n=0 



a * (xq-i/a,xq-i/c;q) c 
(q n+ \aq n+1 /c;qU 



-d n x 



n+j 



(aq n +i+ 1 ,abq n +j+ 1 /c; q)oo {atq^-q). 



« ( g n+l )Cg n+l /Q;g)c 



n+j 



n=0 



(cq n +i+ 1 ,bq n +i+ 1 ; q)^ {ctq n + l ;q) Q 



aqi(aq/c,q;q) c 



(aqi +1 ,abqi +1 /c, atq; q) c 
_ eg 3 (eg/ a, q-^)^ 
(cqi+ 1 ,bqi+ 1 ,ctq;q) 00 

The conclusion of the above calculations is that 

V 



aq^ +1 , abgj +1 /c, atg 

ag/c, 

cq j+1 , bqi +1 , ctq 

cq/a, 



9,9 



Qj(t,a,b,c) 



'■' ' (q;q)j(abq 2 i +2 ,q) 

c 



(4.34) 



(bqj+\c^ +1 ;q) 
(c/a, atq; q) 



-3<P2 



(aqi +1 ,abqi +1 /c; 
H — : s 392 



oo 

;q)oo 



loo 



og J+1 , abqi +1 /c, atq 
aq/c, 

cqi +1 , 6g* +1 , ctg 
cg/a, 



9,9 



9,9 



(a/c,ctg;g) c 

Next we apply (12.5.8) in [15] with A = D/(qat),B = aqi +1 ,C = abqi +1 /c,E = abq 2 ^ 2 and 



q,qct 



let D — > and realize that Qj(t, a, b, c) has the representation 



(4.35) 



Qj(t,a,b,c) 



(q;q)j(aqt,q) oc 
It is clear from ()4.35p that Qj(t, a, b, c) — 
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aq 



abqi +1 jc 



abq 2 ^ 2 



+ 



Next we compute the functions Qj(t, a, b, c). We have 



Qj(t,a,b,c) 



1 



-L(P j (x){-xt;, 



Ai • • • Xj 

(-l)J'g(a) (1 - 9y(ag» +1 ,6^' +1 ,c?»' +1 ,a6^ +1 /q g) c 
a(95 9)i(9, abq 2+2 i,c/a, aq/c; q)^ 

I {-xt; q)ooT> J q wi(x, aq j ,bq j ,cq j )d q x 

J cq 

Vq®(cuj> +1 , bqi +l ,cqi +l ,abqi +1 /c; q) 



a{q; q)j (q, abq 2+2 'j , c/ a, og/ c; g) c 
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where 



9(1 



aq 



aq 



cq 



(xq j /a,xq j /c;q) c 
(x,bx/c, xtq-j;q) x 



E 

n=0 



(-atq n+1 , q n+1 ~ j , aq n+1 ^ /c; q) 



00 „n 



(aq n+1 ,abq n+1 /c; g) c 



— cq 3 



E 



{-ctq 



n+l qn+l-j^ C qn+l-j 



'7a; Q) 



00 „n 



n=0 

a(-atq 3+1 ,q,aq/c; q) a 
{aqi +l ,abqi +l /c] g)^ 
c(-ctg- ?+1 ,g,cg/a;g) 00 



(cqi+ 1 ,bqi+ 1 ;< 



(cg«+ 1 ,6g«+ 1 ; 9 ) 00 

0V+ 1 , a^'+Vc, 
-atq 3+1 aq/c 

cq 3+1 , bq 3+1 , 
-eta- 7 " 1 " 1 , cg/a 



■3<^2 



9,9 



-3^2 



It follows that 

Qj(t,a,b,c) 
(4.36) > 



Pqd) 



(q;q)j(abq 2 + 2 i;q) c 
{bq 3+1 ,cq 3+ \-atq 3+1 ; t 



(c/a; g) c 



-392 



(aqi+\ab^+ 1 /c,-ct^ +1 ;q) 00 , 
H /— 7 v 302 

(a/c; gjoo 



agJ' +1 , abq 3+1 /c, 
-atq 3+1 aq/c 

cq 3+ \ bq 3+1 , 
-ctq 3+l cq/a 



9,9 



To simplify (OHD we apply (12.5.8) in p3] with S = ag J ' +1 ,C = abq 3+1 /c,D = -atq j+1 ,E = 
abq 23+2 and let A —* 00. Alternately we may use the result in Exercise 3.8 of [12]. The 
conclusion is that 



(4.37) Q i (t,o,6, c ) 



? C 2 ) t i(_ atg j+i ;g)c 

(9; g)j 



"2^2 



og J+ , abq 3+1 /c 
abq 23+2 , -atq 3+1 



9, -teg 



j+i 



Finally the transformation (|2.23p gives yet another alternate representation, namely 



(4.38) 



Qj(t,a,b,c) 



q&t 3 

(9; q)j 



(-t;q) 00 20i 



aq 3+1 , cq 3+1 



ibq 23+2 



q,-t 



It is clear from (|P7|) or P~38j) that Qj(t, a, b, c) 



0) 

(<?;<?)j 



+ 



Theorem 4.3. TTte functions Qj(t,a,b,c) and Qj(t,a,b,c) associated with the big q-Jacobi 
functional (|4.30p are defined in (|4.35|> and (|4.38|> (or in (j4.37j) ). Moreover we have the follow- 
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ing addition formula: 



(-gas; 9)00 / qa, qab/c -s/t 
{s;q)oo 3 2 V a ^ 2 , -qas 



(4.39) 



E 

i=o 



g,gci 

{—acstyq^^ +l \aq, bq, cq, abq, abq/c; g) 
(q;q) j (abq,abq 2 ;q) 2 j 



ag J+1 , abqi +1 Ic 



g, get 201 



agJ+\cgJ +1 



g, -s 



Proof. We only need to evaluate T s „Q (t; a, b, c). Clearly T sq Q (t; a, b, c) is 



^ (qa, qab/c; q) n ^ (qa) 

(q,abq*;q) n ^ ^(g;g) fc ^ 



E 

n=0 

00 

E 

n=0 
00 



(qa,qab/c;q) r 
(q,abq 2 ;q) n 



■WE 



k=0 



(q;q)k 



t n + k (-s/t;q) n 



+k 



(qa, qab/c - s/t; q) n 
(q,abq 2 ;q) n 



(qct) n i4>Q(-q n s/t; ;q,qat) 



sr^ qab/c - s/t; q) n (-q n+1 as; , 



n=0 



(q,abq 2 ;q) r 



(aqt; g) c 



The g-binomial theorem reduces the last expression to 

(-qas;q)oo ( qa, qab/c -s/t 

— 7 \ 3<P2 I 7 9 

(qat;q)oo \ abq z , -qas 
and the theorem follows. 



g, get , 



4.4 Second g-Addition Formula 

Define the non-commutative operation S y on f(x) = Y^^Lq a n % n by 



n=0 

where = qxy. Recall [12, p. 28] that the non-commutative binomial theorem reads 

n 

(4.40) ( x + y) n = J2 



k=0 



x k y n - k . 
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Theorem 4.4. The convolution identity (|2,10p is equivalent to the non- commutative addition 
formula 



n=0 



(4.41) S s Q (t; q)=J2*l---K Qn(t] q) Qn(s; q) 
Moreover, we have 

(4.42) (J ,{!.., „ . 

Ai---Aj \{xt;q) c 

Proof. The non-commutative binomial theorem and the definition of gives 



1 W 



fc=0 



Multiply (|2.10p by t m s n / (q; q) m (q; q) n and sum over all m, n > we get 

i m s n v -, t m s n 



(q;q)m(q;q)r 



/ t \ — 7 z~Hj m Hj n Xi . . . Xj. 



This establishes the equivalence of (|4.41[) and (|2.10p . Equations (|4.42p follows from Euler's 



formula. 

Example: We consider the Rogers-Szego polynomials, which are defined by 



□ 



(4.43) 



h n (a;q) = ^ 



fc=0 



They have the g-exponential generating function 

h n (a;q) 



n=0 



(q;q)n (t; q)oo(at; q) c 



Therefore 



S s Q (t; q ) = f2^^(t-s ) 



n=0 



(q;q)n 



(t + s;g)oo(a(t + s);q) c 
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The corresponding orthogonal polynomials are Al-Salam-Carlitz polynomials Un(x; q), which 
have the generating function 



(4.44) 



(t, at; q) 



{xt;q) 

The associated functional is defined by 

1 



x <^n q) 



x — (q;q)n 



(4.45) 

Note that 
(4.46) 



Cf(x) 



(1 - q){q,a,q/a;a) c 



(qx,qx/a;q) 00 f(x)d q x. 



A; 



aq>- l {l-q>) 



It follows from f|4.44|) that the generating function of Q n (x) is 



/I Ai • • • \ n Qn(t; 1 



n=0 



1 {qx,qx/a;q) c 
{q;q) n (I - q)(q,a,q/a;q) 00 J a (xy,xt;q)oo 



y 



(y,ay,q)c 



-d„x 



(ay,q)c 
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(a,t;q) 
(ayt, yq/t, q/ayt;q) 



y, t 

q/a 



Q,Q + 



(y,q)c 



2<Pl 



(l/a,at;q) 
q,q/ayt 
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ay, at 
aq 



q,q 



y, ay 

yq/t 



(t,at,q/t,q/at;q) c 

where the last equality follows from the transformation [12\ (III. 32)]) if we replace all the small 
letters by capital ones and apply the parameter identification: 



Therefore 



^] Ai • • • \ n Qn{t'i 



A = y, B = ay, C = yq/t, Z = q/ayt 
(ayt, yq/t, q/ayt; q) 



n=0 



(q;q) n {t,at,q/t,q/ at;q)oo 



00 ,. I y, ay 

lu " 201 ( yq/t 



q,q5/ayt 



{ayt,yq/t;q)oo 

{t,at,q/t,q/at;q)oo sZl 
{ayt; q)^ 



lim (5q; q) 



00 291 



q/t, q/at 

yq/t 



(t,at;g)oo' 

where we used the transformation [121 (M-3)] and lim (5 ^ 1 -(l — 5) J2^Lo a n& n = ^^n-*oo <h. 
By equating the coefficients of y n we get 



(4.47) 



Qn{t;q) 



1 



(q;q) n (t,at;q) 

c 

Summarizing we get the following addition formula. 
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Theorem 4.5. If st = qts then the following addition formula holds 
1 A, , /•'' 



(4.48) 



(t + s,a(t + s);i 



^(-a)V 2 %a)j 



(g; g)j (t, at; q)^ (q; q)j (s, as; q) c 



4.5 Computing First g-Addition Formulas Using Generating Functions 

Recall that the Al-Salam-Carlitz polynomials {U^ (x; q)} have the generating function (|4.44|) 
and the associated functional is in (I4.45|) . It follows from I14.44j) that the generating function 
of Q n {t;q) is 



y^Aj • • • X n Qn(t; q) 



(y,ay,q)c 



n=0 



(q;q) n (1 - q){q,a,q/a;q) 

_ (ay,q)oo , f y, t 
(a,t;g)oo 2 1 \q/a 



1 (gx,gx/a;g)oo 
(xy,xt;g)oo 9 



+ 



(y;g)c 



>oo , / ay, at 

(1/a, at; g)oo \ 



h q/ayt , 



= {ayt.yq/t.q/ay^q)^ ( y, ay 
(t,at,q/t,q/at;q)oo 2 1 \ yq/t 

where the last equality follows from the transformation [12^ (III. 32)]) if we replace all the small 
letters by capital ones and then take the following substitutions: 



A = y, B = ay, C = yq/t, Z = q/ayt. 



Therefore 



) j \i---\ n Qn(t;q) 



y n (ayt, yq/t, q/ayt; q) c 



n=0 



(q;q) n (t, at, q/t, q/at; q)^ 1™ 

_ {ayt;g)oo 
(i, at; g)^' 

Equating the coefficients of y n in the above identity we get 

t n 1 
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y, ay 
yq/t 



q,q5/ayt 



(4.49) 

Similarly we have 

00 

(4.50) ^Ai---A„g n (t;. 



Qn(t;q) 



(q;q) n (t,at;q) c 



(y,ay;q) c 



n=0 



(q;q) n (1 - q)(q,a,q/a;q) c 



1 (qx,qx/a;q) c 
{xy;q)oo 



-(-xt;q) 00 d q x. 
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The right-hand side of (|4.50p is 

(y,ay;q)c 



{a,q/a;q)o, 
(ay,-t;q) c 
(a;q)oo 



,n=0 



3<P2 



(?;?)« (yq n ;q)oo) 



n=0 



(aa r ' 



-taq n ; i 



(<?;<?)n (ayq n ;q) c 



y,o,o 

-t,q/a 



ay, 0,0 



(y,-at;q)oo 

q,q) + —rr, — n — 392 , , 

(l/o;g)oo V ffa,-o< 



In (III. 34) of [12] replace a, 6, c, d, e by A, 1?, C, D, £?C, respectively then let C — > then let 
A — > oo. The result is the three term relation 



B 
D 



DE 

? '~B~y (E/B;q) 
(B,DE/B;q) 



(£;g)oo M ( 5,0,0 



„ 392 V D,qB/E 



+ - 



392 



5,0,0 



(D^/E-qU V DE/B,qE/B 
We now choose B = y, E = ay, D = —t. Therefore (|4.50p becomes 



q,q 



(4.51) 



Ai • • • X n Qn(t; q) 



n=0 



(q;q)r 



Joo 1^1 



g, —at 



fc=0 



Finally we apply the g-binomial theorem in the form 

(*«).- E 7 rrH A-»)' 

^ (q;q)k(q;q)n-k 

and (|4.46p to obtain 
(4.52) 



Q.( t;g ) = ( tg3;g )°° t^® l9l (0;-^;o,-aV). 



Theorem 4.6. VFe Ziawe t/ie identity 

(-s;<?)oo , / 0,-s/t 

(4.53) 



g, at 



1 



— (-osi) g 1 j i9i(0;-tg ;q,-atq ). 



Proof. We need only to show that the left-hand side in (|4.53p is T Stq Qo(t; q). This can be seen 



22 



as follows. By (|4.43j) and (|4T8|) we have 

00 " t n„fc 



rr- „ , x \ - \ - t n a k (-s/t 



t k a k (-s/t;q) k t n {-sq k /t; q) r 



-'/'x / 0, — s/t 

(t;?)pn V -s 



(959) 

q, at 



□ 



5 Sheffer-type Polynomials 

5.1 Moments of Sheffer-type Polynomials 

For the Sheffer type orthogonal polynomials we can compute Qj(t) by the generating function 



(5.1) ^X 1 ...\ n a n Q n (t)y i =CU xt J2 a nPn(^ , 

n=0 ' V n=0 ' / 

where ay's are some suitably chosen constants. 

The cases of Hermite Laguerre, Meixner, and Charlier polynomials do not lead to interesting 
addition theorems because the addition theorem predicted by Theorem 12.11 follow from the 
binomial theorem. We just indicate the corresponding Qj(t) in the addition formula of each 
family of polynomials. 

Indeed in the case of Hermite polynomials {H n (x)}, 

(5.2) A„ = n/2, H n (x) = 2 n P n {x), 
and 

(5.3) Q n (t) = - exp(t 2 /4). 

n! 

In the case of Laguerre polynomials {L"(x)}, 

(-D n 

(5.4) X n = n(a + n), L«{x) = ^-f-P n ( x ), 
and 

(5.5) Qn ( t ) = -(l-t)-<*- n - 1 . 

n! 
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For Meixner polynomials {M n (x; (3, c)}, we have 

(5-6) X n = ~ 2 1)C , M n (x-p,c)= { j^P n (x), 

and 

1-c Y +n (e*-l) n 



(5-7) Qn(i)=(T 



n! 

In the case of Charlier polynomials {C n (x;a)} 

(5.8) K = an, C n (x; a) = (-a)- n P n (x). 
A calculation gives 

(5.9) QAx) = K - ' exp(e' - 1). 

Ji 

In the case of Meixner- Pollaczek polynomials P„ (x; </>), |18j . 

ra(re + 2A-l) X/ (2sin0)" . , 

4sm <p n\ 

One can see that the Qj's are given by 

Note that for the orthogonal polynomials of Sheffer type all the Qo W s have been given in [29|. 

5.2 Sheffer-type Polynomials as Moments 

For the Hermite polynomials, we have 

(5.12) ew (2xt-t 2 ) = jr^^t n . 

n=0 

Let 

Qn{t) = -e 2xt - 1 = - + 2(n + ry + ■ ■ ■ . 

n! n! (n + lj! 

From exp (2x(t + s) — (t + s) 2 ) = exp (2xt — i 2 + 2xs — s 2 — 2is) we derive the addition for- 
mula. 
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Theorem 5.1. The functions {Qj(t)} have the addition formula 

oo 

(5.13) Q (t + s) = ^2 n\{-2) n Q n {t)Q n {s). 

n=0 

It follows that A n = — 2n and 

. )H n (x). 

Remark. Radoux [22] proved (|5.13p by computing the corresponding Stieltjes tableau using 
induction. 

For the Laguerre polynomials Ln {%) we have [18\ p. 48]: 

(5.15) e* P1(-; a + 1; -xt) = £ ^. 

In (|3.8p letting £> = a + /? + 1 and substituting i and s by t/a and s/a, respectively and, then 
let a — > oo we get 

5 + 1; -2(t + = £ VgwifiT o^H B + 2n + 1; -2t) 

(5.16) ^ (^) 2 n( J D + l)n 

x F 1 (-;B + 2n + l;-2s). 

Let 

t n 

(5.17) Qn(t;a) = — e* Fi(-; a + 2n + 1; -2arf), n > 0. 

n! 

Then we have the following addition formula. 

Theorem 5.2. The functions {Qj(t;a)} have the addition formula 

™ nUcv) (-4r 2 ) n 

(5.18) Q (t + S ; a) = V /\ 7 ^ Qn(*i a)Qn(*; a). 

~ (a)2n(« + l)n 

As an immediate consequence we have A n = ( a+ 2n-i)™a+2n-2) (L+n) anc ^ 



(5.19) H 



(a + 2i + l) n 
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For the Meixner polynomials M n (x; (3, c) we have 

1 - c 



(5.20) 



e l 1 F 1 







t j = yv M n (x;0,c) tn 



n=0 



rv. 



In (|3.8p substituting a + 1, a + /3 + 2, i and s by — x, /3, (c— l)i/2c and (c— l)s/2c, respectively, 
we obtain 



(5.21) 



/3 



1 — x 



n=0 



n!(-x) n (/3 + x) n (/3 - l) n (I - c 



(/?-l) 2n (/?), 



2n 



(t s y 



n — x 
l3 + 2n 



1 - c 



t) 1F1 



n — x 
/3 + 2n 



1 -c 



1 -c 



Therefore define 

Q„(x;/3,c) = ^ e t 1 F x ( 
we have the following addition formula. 

Theorem 5.3. The functions {Qj(t; (3, c)} /iaue i/ie addition formula 



(5.22) Q Q (i + S ;/3,c) = J] 

n=0 

In the same way we derive 



n\(-x) n (/3 + x) n (/3 - l) n fl-c 
(/3-l) 2n (/3) n 



2/i 



Q n (t;/3,c)Q n (s;l3,c). 



A, 



n(-x + n - l)(/3 + x + n- + n - 2) /"l - c 
(/? + 2n - 2)(/7 + 2ra - 3)(/3 + ra - 1 



and 
(5.23) 



if, 



,i+n 



i + n 



M n (x-z;/? + 2z,c). 



The Meixner- Pollaczek polynomials P^{x;(f)) have the generating function 



(5.24) 



e*iFi 



A + ix 
2A 



1)t)= f |W) ( , 



n=0 



(2A)„e in * 



In ((MJ) substituting a + 1, /3 + 1, t and s by A - 1 + ix, A - 1 - ix, (1 - e _2i<?i )i/2 and 
(1 — e~ 2lr f')s/2, respectively, and letting 



(5.25) 



QW(x;<f>) = -e\F 1 



in 



\ + ix + n 
2A + 2n 



-2id 



1)* , 



we obtain the following addition formula corresponding to Meixner-Pollaczek polynomials. 
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Theorem 5.4. We have the addition formula 

(5.26) Q?>( ( + ,*) = f; "KA + i,MA -i,)„(2A-l)„ 4 „ Q!iX , (t; ^ 

6 g-Ultraspherical and Askey— Wilson Polynomials 

The continuous q-ultraspherical polynomials have the weight function 

1 (e™,e-™) (/3 2 ,g)oo 1 
M^P) = 77- /^ 2i8 o -2iB\ (a a \ H 9 ' x = cos 6 6 (—1, 1), 

and have the property 



(«)„ " J 4(l-/3«J-i)(l-/3^) 
In view of (jl.lOp and (|1.8p we see that 



2 °° 



(6.1) e^ = -^(n + l)/ n+1 ( y )t/ n ( 



^ n=0 



Now the special case 7 = q of the connection relation (jl.9p gives the following expansion: 

k (q/P;qMg;q)n-ki-Pq n - 
{q;q)k(q(3;q)n-k 1-/3 



(6.2) [/„(*) = ^ T^w^V: . ——5 C„_ 2fc (x;/3 I g). 



fc=0 

Using the above expansion and the orthogonality of the g-ultraspherical polynomials we have 
2 1 00 ri 

Qj(t) = 77 r J2 / w ( x )( n + i)4+i (i)^(^) Pi(x)dx 

{q)j {q;q)k{qP;q)j+k l-P 

To denote the explicit dependence on g and /? we set 

^on n / + * „n 2i+1 , 1U m /3 fe (g//3;g)fc(g j+1 ;g)fe 

(6.3) Qj{t;P,q) = — — 2^0 + 2fc + l)/ j+2 fc+i W" 



t ^ v ;.™w (</;<7)*(pV +1 ;<z)* 



Thus we proved that 



27 



Theorem 6.1. The functions {Qj(t; f3,q)} have the addition formula 
(6.4) 



Q (s + 1; fi, q) = f; }j^ iP ^ q h Qn (* = > <l « <■ </)■ 



n=0 



4»(j3',q) n (ql3;q) t 



The special case /3 — > is worth recording. Indeed if 

nn+l 00 



(6.5) 



F n (t; q) := — Y;( n + 2k + l)In+2k+l(t){-l) n q 



k=0 



n + k 
k 



then we have established the curious result 

(6-6) F (s + t; q) = ^^F n (s; q)F n (t; q). 



n=0 



Another interesting case is to let f3 = q v then let q —>■ 1. This should reduce (|6.3p to (13. 1J) 
since lim 9 _>i C n (x; q^lg) = C^(x). Surprisingly the q — ► 1 limit of (|3.ip . after setting /3 = is 



Q j (t) = 2^I i+2 ,(t) 



fc=0 



(-v)k(j)k 3 + 2k 
k\{j + l + v) k j 



Equating the above limit and the Qj as in (J673J) leads to the following known identity involving 
Bessel functions 



(6.7) 



n=0 



n! r(/i + 1 — /x — n)r(z/ + n + 1) 



see (7.15.2)]. It is also worth mentioning that (|6.7p is equivalent to a theorem of Bailey 
evaluating the sum of a well-poised 4-F3 with argument —1, [6j (4.5.4)] . 

Next we consider the Askey-Wilson polynomials whose weight function is 



W(x;ai,a2,az,a 4: \q) 



(e 2id ,e- 2id ;, 



1 



(6i 



(g;g)oorii<i<fc<4( a i a fc;9)c 



X = cos ( 



2tt (aia 2 a 3 a4; g)^ 

The Askey-Wilson polynomials have the basic hyper geometric function representation 

p n (x;ai,a 2 ,a 3 ,a 4 | q) = a^ n (aia 2 , aia 3 , aia 4 ; <?)„ 

, { q~ n , a^aza/tq™' 1 , aie 10 , aie~ l6 

X 4</>3 ( ______ -> 9 



(6.9) 



ai<22, 0103, 0104 
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One very special case of their connection coefficients formula is [21 (6.4)— (6.5)] 

n 

(6.10) p n (x;a,a 2 ,a 3 ,a4\q) = 

(x; a, a 2 , a 3 , 04 | q) 



k=0 



where 



a"- fc (g; 9)n(«a 2 a 3 a4g"- 1 ; g) k (a/a; g) n - k ^ , k , 

(6.11) c kn = -, 7— ; — r— ; gfc — \ [ajamq ;q) n - k . 

(q,aa 2 a 3 a4q k: i ;q) k (q,aa 2 a 3 a4q Zk ;q) n - k 2< <^ <4 

Moreover 

^(^) = ( g n+2. g ) Pn(g;Vg!g>~Vg)~g|g)j Pn(x;a,b,c,d\q) = 2 n (abcdq n ' 1 ;q) n P n (x) 



(1 - g n )(l - aia 2 a 3 a 4 g n 2 ) ]li<j<fc<4( 1 ~ a j a kq n ^ 



4(aia 2 a 3 a4g 2n " 3 , aia 2 a 3 a4g 2n - 2 ; q) 2 
Applying ()2. 18|) and the plane wave expansion (|6.ip we find that the Q m 's are given 

Qm(t) = t— : ^C(e xt P m ) = - - 1 - V(n + l)/ n+ i(t)£(f7 n (z)P m ) 

A1A2 • • • A TO A1A2 • • • A m ^— * 

71=0 



= (2/t) £ (n + i)/n + i(t) ( ;C 2 ! ; ? m 2n W 

n=m ^ ' 

fou\^^f , 1 \ r ^ ( Q g m+1 ; g)m a n ~ m (g; g)n(g n+2 ; g) m (g/a; g) n - 

(2/t) ^2 (n+ l)in + i(t) 2- 



(9 n+2 ;g)n (g,ag m + 1 ;g) m (g,ag 2m + 2 ;g) n „ m 
x(-9 m+1 ,? m+3/2 ,-g m+3/2 ;g)n- m - 
After some simplification we arrive at 

nm+l 00 

Qm(t) = — ^ + m + !) 

(6.12) n=0 

^ (q m+1 ,q/a,-q m+1 ;q) n (g 2m+3 ;q 2 ) n 

(q,aq>™+\q n +™+ 2 ;q) n in+ra+lW ' 
Theorem 6.2. TTie functions {Q m (x)} defined in (|6.12|) satisfy the addition theorem 

00 1 2 2 \ 

(6-13) Q (S + t) = J2 ,n7 ' 2 ? Qn(t)Qn(s). 

^^ n [aq,aq l ;q) 2n 
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7 Ultraspherical Polynomials as Moments 



One of the generating functions reads 
(7.1) 



00 cfth-A 



(2u% 



(x 2 - l)V 



v + 



1 ) ' 



n=0 ^ '" v ■ 2 



Let coscp = -1, then w = 2 + Zand C£(-l) = (-l) ni ^. It follows from (TT5]) that 



0-^1 



(7.2) 



u 


+ 1/2 




n 




0-^1 



-(z + Z) ! 



E 

n=0 
~2 



(w + i/-l/2)(-l) w (2t/-l) T 
n!(i/ + l/2) n (i/ - l/2) n+ i 



z/ + l/2 + n' 4 7° l \u + l/2 + n' 4 
Therefore, let z = ty/l — x 2 and Z = ty/l — x 2 we obtain 



(7.3) 
where 



Q(f + .) _ f <"r:ffl ( : 1) : ( *- 1M '-^W (.) 



71=0 



Qn(t) = - 1 e tx F 1 



{u + l/2) n (u-l/2) n+1 



(x 2 - l)t 2 



t 



t 



n+1 



n\ V v + 1/2 + n ' 4 

Extracting the coefficients of t m s n in (|7.3p we get 



(m + n)!C£ +n (s) _ ~ (fe + !/ - l/2)(-l) fe (2z, - l) fe (1 - x 2 ) 



(7.4) 



(2i/) m+n m!n! 



E 

fc=0 



+ (n+ l)x- 
n! (n + 1)! 



2\fc 



+ 



fc!(i/ + l/2) fc (i/ - l/2) fc+1 4 fe 



Using the relation 

we derive 
(7.5) 



(2^) m -fc(2^) n _fc 

lim a-"/ 2 C« +1 / 2 (x/v^) = ^4^, 



/?,! 



m!n! 



mAn 

E 

fc=0 



(-2) fc g m _ fc (a) g ra -fe(a) 
fc! (m — fc)! (n — A;)! 



An immediate consequence of (|7.3p is the following formula for the Hankel determinant eval- 
uation \28\ Corollary 3]. 
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Corollary 7.1. We have 



and more generally, for n > 0, i/ie entries of the Stieltjes tableau (|2.9|) are 

(7.7) ff M+n = £ (x 1} 1 



i\k\(n-2k)l(u + l/2 + i) k 4 k ' 



8 A variation of the Stieltjes-Rogers addition formula 

Let {P n (x)} satisfy (|2.ip with moment sequence {[J, n } and P n {x) = a~~ n P n (ax + b) (a ^ 0). 
Then it is well-known j5j p. 25] that {P n (x)} is an OPS with respect to the moments given by 

(8-1) 0n = a- n E(?)(-&) n- V*, 

k=0 ^ ' 

and satisfy 

(8.2) P n+1 (x) = (x - b -^) P n (x) - %P„-i(x). 

\ a J a z 

Let Qo(t) = J2^L Mnfrr- Then it is easy to see that 

(8.3) Q (t) = V fin- = e- bt / a Q (t/a). 

£ ' 77 



n=0 



The following variation of the Stieltjes-Rogers addition formula (12. 12ft is sometimes very useful. 
Theorem 8.1. The addition formula for the moment sequence (|8.ip is 

oo 

(8.4) Q (s + t) = ^ Ai . . . X n a- 2n Q n (s)Q n (t), 

n=0 

where Q n (t) = e~ bt ^ a Q n (t/a). The corresponding entries in (|2.9p are 

n+j 



i.5) ^• + n=E( n r J )(- 6 )* o " a " ifl i, 

fc=0 ^ ' 



31 



In particular we have H n ,n = (—b/a) n H njn , i.e., 

(b\ n 
— j det(ni + j)o<i,j< n - 

For example, let fi n = (a + l) n be the nth-moment of Laguerre polynomials {L®(x)} (see 
U53D). Ha = b = l/x then 

(8-7) /2 n = ^(-l)"- fe (^x fe (a + l) fe 

fc=0 ^ ' 

is a weighted derangement number. The corresponding addition formula reads 

oo 

(8.8) Q (s + t) = J2 n\(a + l)nX 2n Q n (s)Q n (t), 

n=0 

The a = case of the above formula was derived by Radoux [22J using induction. 
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